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ABSTRACT 


The analysis of vibration generation and transmission in shipboard 
machinery installations has defied rigorous mathematical description. 
After certain simplifying assumptions have been made, the vibration response 
of the installation can be predicted analytically. Theeobject of this thesis 
is to demonstrate the method by which the simplified system can be analyzed 
to give vibration trends and characteristics. The method incorporates the 
four pole parameter technique of analysis and a systematic study of the 
effects of installation modifications. 


The four pole parameters of a system are those expressions which 
inter~relate two basic input and output quantities. For vibrating mechanical 
systems these quantities are force and velocity. The four pole parameters 
must be derived or measured for the mechanical elements which transmit the 
force and velocity through the installation. In a shipboard installation 
these elements are the machine, mounting, foundation, hull, and water. 

After the four pole parameters are obtained for these five elements, 
matrix notation and combination permit a transmission study of the entire 
five element system. The computational assistance of a digital computer 
is an indispensable aid in this portion of the study. 





Numerous conclusions and recommendations were obtained from the 
analysis and from the results of the systematic study. The most important 
of these were? 


(1) For the purpose of vibration analysis a complex machinery 
installation can be represented by a relatively simplified 
system. This representation and the use of the four pole 
parameter technique are sufficient to establish trends and 
characteristics. 


(2) The stiffness of the resilient mounting has a dominant effect 
on the vibration transmission at all frequencies. 


tal 





(3) A Portsmouth mounting with decreased stiffness and increased 
damping and synthetic rubber damping applied to the foundation 
were the most effective installation modifications in providing 
overall transmission reduction. 


(4) As new ideas for reducing vibration transmission are intro- 
duced, they should be investigated by means of the four pole 
parameter technique. 


(5) A model of the approximate system of this thesis should be 
constructed to permit experimental verification of the four 
pole parameter technique. 


Thesis Supervisor: John R. Baylis 


Title: Associate Professor of Naval 
Ingineering 
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CHAPTER I 
INTRODUCTION 

1.1 Importance of Vibration Transmission 

The problem of noise control is one of major importance to naval 
engineers. A ship with a high noise level operates under the joint 
handicaps of being easily detectable by the enemy's listening devices and 
of reducing the effectiveness of its own sonar. Noise control implies 
the control of vibrations since acoustic or noise energy in a medium is 
synonymous with vibration transmission in that medium. A thorough under- 
standing of the vibration transmission characteristics of shipboard 
installations can provide the solution of the noise control problem. 
1.2 Current Developments 

Spurred by the Navy's need for quieter submarines, many facilities 
and laboratories have been devoting time and effort towards finding a 
solution of the noise problem. Many studies are currently in progress 
seeking mechanical means of reducing vibrations in machinery installations. 
Much of this is being financed by Navy contracts; thus, most of the major 
new developments have shipboard applications. 

The development receiving the most application at present in the 
Navy is the rubber-type resilient mount. These Navy mountings have been 
developed by the Engineering Experiment Station, Annapolis, Maryland and 
the Portsmouth Naval Shipyard, Portsmouth, New Hampshire (15). The 
effectiveness of resilient mountings has been studied by Mr. A. O. Sykes 
of the David Taylor Model Basin (13, 14). Studies have also been made of 
the properties of other rubber-like materials and of their effectiveness 
when used as mountings (17). 

It is becoming more apparent to interested engineers that the 


other elements in the vibration transmission path from source to water must 


(1) 





(2) 
be investigated. Work has been done in determining vibration character- 
istics in machinery, the most common vibration source (12). Machinery 
foundation vibrations and techniques of damping them are being experimen- 
tally evaluated by the Sound and Vibration Section at the Electric Boat 
Company (20). ‘This work has emphasized the importance of the foundation 
in the vibration transmission path. 

The final item in the vibration transmission path is the ship's 
hull. Measurements of the vibration response of submarine hulls have 
been made. This empirical data gives a ane means of predicting the 
interacting effects of foundation, hull, and water (11). 

It has been recognized that the time to apply noise control 
measures aboard ship is in the preliminary design and construction stages. 
This concept has been applied to the SSN 593, Thresher, to be commissioned 
this summer. A consistent set of specifications, setting allowable noise 
and vibration levels, have been developed for this ship. 

1.3 The System Under Study 

There are several means by which noise or vibration aboard ship 
can be transmitted from its source to the water. If an operating machine 
is considered as a typical noise source, its vibration can reach the hull 
and water by several paths. These are airborne transmission, supporting 
structure transmission, and transmission through connected piping. The 
vibration paths are best illustrated by the component block diagram in 


Figure I. 





(3) 
Figure, I 


System Block Diagram 


Machine Mounting Foundation 





The initial step in formulating an analytic approach to this 
transmission problem is to assume that each of these paths can be treated 
separately. This thesis is limited to an analytic study of the trans- 
mission through the supporting structure. The path under study, as 
indicated in Figure I, consists of machine, mounting, foundation, hull, 
and water. 

1.4 The Analytic Approach 

The exact describing equations of the transmission phenomena 
through the above mentioned mechanical elements would be an ideal analytic 
tool for the investigation of this problem. The actual system, however, 
defies any clear-cut description, and to obtain the describing equations 
is a Herculean task. The system is non-linear and possesses many degrees 
of freedom. The vibration of most machines is of a random nature. In 
addition, these machines are mounted on complex structures. 

To circumvent the above complexities, simplifying assumptions 
must be made to obtain a workable mathematical model. These assumptions 


are completely enumerated in Chapter II of this thesis. The simplified 





(4) 
model is not intended to give the complete picture of vibration trans- 
mission in the actual installation. 
The object of this thesis is to make a systematic study of the 
analytically-described, simplified system. It is hoped that from this 
simplified study vibration trends and preferred design practice for the 


complex machinery system will become available. 
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CHAPTER IIL 
PROCEDURES 
eel The Preitamimery Experiment 

The authors first became interested in an analytic approach to 
vibration transmission in shipboard installations while on temporary duty 
at the Admiralty Engineering Laboratory, West Drayton, England. A series 
of experiments were being conducted there to determine isolation charac- 
teristics of Admiralty resilient mountings. Results had been obtained 
from one of these basic experiments, and the aGeHexe investigated various 
approaches to the mathematical verification of this data. 

The experimental results studied were obtained by driving a beam 
Supported by an Admiralty mounting with a typical vibration generator. 
This experiment is described in Appendix A-l. 

Most previous experiments investigating the characteristics of 
mountings treated the supported member as a lumped mass. This treatment 
facilitated use of a theory involving electrical analogies (14). In this 
case the beam was used as the supported member because of its non-rigidity. 

A recent paper by R. D. Cavanaugh and J. E. Ruzicka provided a 
method by which a vibrating beam-spring system could be analyzed (4). 
This method was the four pole parameter technique. The four pole param- 
eter method had been widely used to analyze mechanical systems; therefore, 
it was selected as the technique to provide the mathematical verification 
of the beam-mounting experiment. 

2.2 ine oun rove Parameter Technique 
eee Seocnera. 

There are many systems in nature which have two basic quantities 
as input and two basic quantities as output. An electric circuit element 


is one of these systems. In this case voltage and current are both input 





(6) 
and output. All these systems are called "four poles" since four quantities 
become inter-related by some action within the system. 

Klectrical engineers have devised a technique of analysis for 
eircuit elements called the four pole parameter technique. Each circuit 
element is considered independently, and complex circuitry is developed 
by connecting the various elements. Recently the four pole parameter 
technique has been applied to mechanical vibration problems where the 
mechanical element considered has a force and velocity as both input and 
output. 

A typical mechanical system can be represented in block diagram 


form as illustrated in Figure II. 


Figure Il 







Mechanical 


system 





FL and Vv, are the input force and velocity, while Fy and V5 hss qrlahs 


output force and velocity. This system can be described in general by 


the following equations; 


z 1 
Fs C44 PF, “ E15 Vy ( ) 

5 2 
ty Ga ay ap os (2) 


21? and € are the four pole parameters of the mechanical 
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system. 


(7) 


Using matrix notation, equations (1) and (2) become: 


a ane - “12 By 
= (3) 
a Sn coe Vo . 


——- 


The determination of the four pole parameters comprises the ma jorportion 
of the technique. Various rules must also be developed to permit combina- 
tion of the matrices to represent mechanical elements in various multiple 
configurations. 
Pee wc Development of the Four Pole Parameters 

Basically the solution for the four pole parameters of the mechanical 
system can be reduced to the following procedure: 

1. Investigate the system and obtain its performance equations. 

2. solve the equations imposing the proper boundary conditions. 

If non-linearities exist in the system it may be necessary to 

apply linearization techniques at this step of the procedure. 

3. Reduce the solutions to the form of equations (1) and (2). The 

four pole parameters can then be written by direct comparison: 

This procedure is applied to several basic mechanical elements in 
Appendix Ame. 

The above procedure can be specifically applied to the mechanical 
systems found in the analysis of vibration transmission. Equations (1) 
and (2) can be solved for the four pole parameters by imposing the condition 
that the output, which is station 2 in Figure I],is either blocked or 
free. The blocked condition implied the velocity at station 2 is zero; 


the free condition implies the force at station 2 is zero. 


with station 2 blocked (4) 


cal he 
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F 
E _ with station 2 free (5) 
oe V5 
“a 
E —— with station 2 blocked (6) 
Ol Fy, 
vs 
Bp ia with station 2 free Cp) 
Co Vy 


The following definitions can be applied to equations (4) through 
(7) to simplify the notation. 


F 


(T 2 = Perce Imsasicslen == with station 2 blocked 
F 
al 
ata Vo 
ed) = Velocity Transmission = >= with station 2 free 
, ; 
Bi Vo 
(M, .) = Transfer Mobility =—— with station 2 free 
2 
b ia 
ul) = Transfer Mobility === with station 2 blocked 
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Equations (4) through (7) become: 


e414 = hy (4a) 
E15 = i) (5a) 
“51 (M4) (6a) 
S29 = (R (Ta) 


Equation (3) can now be written: 
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This development provides a convenient notation by which any 


system may be analyzed. 
ee. Combination of Four Pole Matrices 
Multiple mechanical elements for which the four pole matrices 


have been derived can be connected in series. This is illustrated in 


fieure Lil. 


Bicure Lit 





If this is the case, the overall system of n subsystems Sy through Sn 


can be represented by a single four pole matrix [S_]. 


The single four pole matrix [s.] is defined: 


[SJ = [S, ] x [S,] x iS] x cocerereeeelS) -~1]x [s | (9) 


Multiple elements can also be connected in parallel as illustrated 





ineragure IV. 





These n parallel subsystems can be represented by a single four pole 


matrix [S ]. 
Pp 


aT Bo 
= om 
ni Yo 


The derivation of the matrix is. is more complicated for the parallel 

than the series case. The solution for the four poles of is, is included 

in Appendix A-3. 

2.2.4 Application of Four Pole Parameter Technique to Preliminary Experiment 
From the preliminary experiment mentioned in section 2.1 and fully 


described in Appendix A-l, experimental results were available to test the 
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validity of the four pole parameter technique. Using this technique the 
isolation, which is the ratio of input force to output force for the systen, 
was theoretically obtained as a function of driving frequency. The experi- 
mental results closely agreed with the values of isolation determined by the 
four pole parameter technique. This demonstrated the effectiveness of the 
technique. 
2.3.1 Description of the General System 

As mentioned in Chapter I, the overall system which is the subject 
of this thesis consists of the mechnaical elements in the structural path 
of vibration transmission from the machine source to the water. The 
block diagram of this system is Figure V. 


Figure V 


Machine Mount ing Foundation Hull Water 
1 2 3 4 15 6 


2.3.2 The Analytic Study 
The initial step in any analytic study is the enumeration of the 
simplifying assumptions. In this thesis these assumptions are: 
1. The transmission paths in Figure I can be considered separately. 
2. The principal contributors to structural vibration transmission 


are the force and velocity normal to the hull. All others are of 


(a2) 

lesser importance. (This assumption has been corroborated by other 

investigators (17).) 

3. The mechanical elements are elastic and linear. 

4. The vibrations in the system are periodic. 

Other assumptions have been made when analyzing the individual elements. 
These will be discussed in deriving the four poles of each mechanical 
element. 

It is recognized that the above assumptions will not permit an 
exact solution of the vibration transmission characteristics of an actual 
machinery installation. As stated in Chapter I, however, the purpose of 
this thesis is not to accurately predict vibration transmission for any 
particular mechanical system. Rather it is to provide a systematic 
study of a hypothetical installation such that trends caused by parameter 
variations can be predicted. A knowledge of these trends could indicate 
design procedures to reduce vibration transmission. 

With the use of the simplifying assumptions and the equations of 
dynamic analysis, the four pole parameter technique can be applied to 
the general system described in Section 2.3.1. The dynamic analysis 
relates input force and velocity to output force and velocity for each 
element. The principal task is the determination of the four poles which 
appear in these equations. 

2.4.1 The Machine 

The problem of mathematically representing the vibration charac- 
teristics of an operating machine is a very difficult task. A study to 
evaluate the internal impedance of a typical machine has been made by the 
A. D. Little Co. (12). In this report an ideal machine is considered to 


possess many resonances and a complex response spectrum. Even though in 
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some cases the exciting internal forces are random, many machines rotate 
and thus exhibit a form of periodicity in their force generation. For 
this reason a harmonic force generator is considered as the vibration 
source in the A. D. Little report. In their mathematical model a harmonic 
force generator drives a multi-mass-spring system in which the relative 
spring stiffnesses and mass sizes determine the resonances. Some investi- 
gators have concluded that this assumed lumped parameter behavior is valid 
at frequencies of only 150 eps or less (17). 

Recognizing the limitations of the lumped parameter treatment, 
this thesis postulates that the machine can be represented by a simple 
beam driven by a harmonic force generator. A simple beam by virtue of 
its structural properties possesses distributed parameters and eliminates 
the limitations imposed by the lumped parameter analysis. 

The four pole parameters for the beam are developed in Appendix 
Fah). 

2.4,2 The Mounting 

Of the mechanical systems in the structural transmission path, 
the mounting is the one best approximated by theory. Historically the 
installation of rubber-type resilient mountings was the first remedial 
step taken to reduce vibrations. Consequently, more mounting theories 
have been developed and tested than theories for the other mechanical 
elements. 

The theories used to develop the four pole parameters of the mounting 
in the general system are those proposed by D. V. Wright and A. C. Hagg 
Given), &. ©. Sykes (13, 14), and/J. C. Snowden 167 

The first mounting considered in the general system is the standard 
Portsmouth mount, which is found in many shipboard installations today. 


The describing four poles for the mount are derived in Appendix A-4.2 
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according to the theory developed by Wright and Hagg. 

The second mounting considered in the general system is one 
constructed of a resilient rubber-type material similar to that bearing 
the trade name of Thiokol R. D. The four pole parameters for this mount 
are developed following the concepts of Sykes and Snowden in Appendix 
A-4.2. 

2.4.3 The Foundation 

In most practical installations, machinery foundation structures 
are combinations of simple beams. This fact naturally suggests the use 
of beam vibration analysis to describe the foundation element. This 
analysis has been used by a number of investigators as a basic step in the 
description of the dynamics of foundation vibration (17, 20). 

For the purpose of this investigation the foundation was simplified 
to be a beam resting on two columns. The describing four poles for this 
mechanical system are derived in Appendix A-4.3. 

2.4.4 The Hull 

The hull of a ship is a complex structure. Unlike other elements 
in the structural vibration transmission path, it cannot be represented by 
any single mechanical element or simple combination of elements. Experi~ 
mental results are necessary to describe the response of typical hull 
sections (7, 11, 17, 20). 

A mathematical model devised by D. V. Wright has provided a 
reasonable approximation of the vibration characteristics of a ship's hull. 
This model is used to obtain the four poles of the hull section in Appendix 
A-4 4. 

2.4.5 The Water 
The water is considered only in its effect as an attenuation 


element. Any entrained mass effect is included as part of the hull. The 
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water presents an effective specific resistance to the propagation of 
logitudinal pressure waves. This specific resistance is dependent on the 
density of sea water and the velocity of sound in that medium. 

The describing four poles are derived in Appendix A-4.5. 
Bap. SOlcion of the Problem 
0 eae General 

The derivation of an expression for vibration transmission 
through the mechanical elements of the general system is the key to 
achieving the purpose of this thesis. The solution of this problem 
involves the combination of the five basic subsystems described in Section 
2.4 into a single complex system. The four pole parameters of the combined 
system are then obtained according to the equations of combination explained 
in Section 2.2.3. The vibration transmission is defined as a force ratio. 
The force transmitted into the water is considered the significant noise 
producing force in the system. Therefore, the vibration transmission is 
defined as the ratio of the force into the water to the input force. 

Since many of the four poles in the system are functions of fre- 
quency, the final ratio will be a function of frequency. Evaluation of ° 
this ratio at sufficient frequencies will give the desired transmission 
versus frequency data. To conform to convention transmission is converted 
to decibels. 

The combination of the general system fourppoles and a sample 
ealculation at one value of frequency are included in Appendix A-6. 
PPS ys Computer Application 

This thesis could not have been completed without the computational 
assistance of the IBM 709 digital computer. The solution for vibration 
transmission at each frequency required a calculation of length comparable 


to that of Appendix A~6. This study required thousands of such calcula- 
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tions. The matrix notation of the four pole parameter technique was 
easily adapted to digital computer methods. 

A computer program was written to solve the combination of five 
four pole matrices. Provisions were made to change any of the mechanical 
elements in the system of Figure V. A complete sequence of this computer 
program logic is explained in Appendix A=-7. A most convenient feature of 
the program was the use of the cathode ray tube oscilloscope to graphically 
present the results. The scope photograph eliminated the necessity of 
manual plotting. 

226 Systematic Study to Determine Trends 

With the selection of the general system of Figure V and the 
mechanical elements in Section 2.4 to represent this system, the problem to 
be considered was established. The four pole parameter technique described 
in Section 2.2 provided an analytical means to study the problem. The 
final phase of the study, however, was the most important. It was to 
devise a consistent, systematic method of changing the general system 
such that many different cases might be evaluated and compared. This 
comparison and evaluation formed the basis for determining the results, 
conelusions, and recommendations of this thesis. 

A detailed description of the systematic changes to the general 


system is given in Appendix A-=5. 
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CHAPTER IIT 
RESULTS 
3.1 The Photo Plots 
The results of the analytic study of vibration transmission in 
shipboard installations are presented in the form of forty-one photo plots. 
These transmission versus frequency plots are included in Appendix B with 
one plot for each of the cases mentioned in Appendix A-5.7. 
3-2 Statements of Results 
The following statements of results are deduced from a study of 
the photo plots. To preserve the systematic procedure of the study, each 
case is examined separately. For detailed descriptions of the systems 
studied, refer to Appendix A~5. 
Case l. This was the basic case involving System lL. It was 
selected as the basis for comparison since the elements comprising 
system 1 were typical of those presently found in shipboard 
installations. 
Case 2. The System was the same as Case 1, except for the 
following: 
a. Mount damping increased. 
The first transmission peak was lowered. 
Case 3. The System was the same as Case 1, except for the following: 
a. Mount stiffness increased. | | ; a. 
The transmission over the frequency range studied was increased. 
Case 4. The System was the same as Case 1, except for the following: 
ae Mount stiffness decreased. 
The transmission over the frequency range studied was lowered by 
approximately 10 db. 


Case 5. The System was the same as Case 1, except for the 
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following: 
ae Thiokol R. D. mount replaced Portsmouth mount. 
The first transmission peak was lowered while the second peak 
was removed. 
Case 6. ‘The System was the same as Case 1, except for the following: 
ae Hull damping was doubled. 
There was no change from Case l. 
Case 7. The System was the same as Case 1, except for the following: 
a. Mount damping increased. 
be Mount stiffness decreased. 
The results were the same as Case }. 
Case 8. The System was the same as Case 1, except for the following: 
ae Thiokol R. D. mount replaced Portsmouth mount. 
b.- Hull damping was doubled. 
The results were the same as Case 5. 
Case 9. The System was the same as Case 1, except for the following: 
a. Mount damping increased. 
be Mount stiffness decreased. 
ce. Hull damping increased twenty times. 
The transmission was two or three decibels below that of Case 7 at 
the lower ee 
Case 10. The System was the same as Case 1,, except Tor the 
following: 
a. Viscous foundation damping was added. 
The initial transmission slope was changed. The transmission was 
lowered at frequencies greater than 1700 cps. 
Case ll. The System was the same as Case.1, except for the 


following: 
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a. Thiokol R. D. and Portsmouth mounts in parallel. 
The transmission was increased at frequencies greater than 1000 cps. 
Case 12. The System was the same as Case 1, except for the 
Pollowing: 

ae Water resistance varied with frequency. 
The transmission was lowered at frequencies greater than 1700 cps. 
Case 13. The System was the same as Case 1, except for the 
following : 

ao Viscous foundation damping was added. 

b. Mount stiffness decreased. 
The first transmission peak was shifted from 21 eps to ll cps. 
Transmission was decreased over the entire frequency range by 
at least 10 qb. 
Case 14. The System was the same as Case 1, except for the 
following: 

ae Thiokol R. D. foundation damping was added. 
Transmission was lowered appreciably at frequencies above 800 eps. 
Case 15. The System was the same as Case 1, except for the 
fowLowing : 

ae Thiokol R. D. and Portsmouth mounts in parallel. 

b. Portsmouth mount stiffness decreased. 
The results were the same as Case 4 to 500 eps. Above 500 eps, 
the results were the same as Case ll. 
Case 16. The System was the same as Case 1, e: cept for the 
LO} owed: 

a. Huil damping increased, maximum. 
The transmission was reduced over the entire frequency range. 


Case 17. The System was the same as Case 1, except for the 
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following : 

a. Thiokol R. D. mount replaced Portsmouth mount. 

b. Thiokol R. D. foundation damping was added. 
The first and second transmission peaks were removed. The 
transmission was reduced at frequencies above 800 eps. 
Case 18. The System was the same as Case 1, except for the 
following: 

ae Thiokol R. D. foundation damping was added. 

b. Mount stiffness decreased. 
The vibration transmission was reduced 10 db. in the 10 to 800 eps 
range. Above 800 eps, the transmission was further decreased. 
Case 19. The System was the same as Case 1, except for the 
following: 

ae Thiokol R. D. foundation damping was added. 

be. Mount stiffness decreased. 

ce Mount damping increased. 
The results were the same as Case 18, except the first peak was 
removed » 
Case 20. This was the basic case involving System 2. Foundation 
mass increase was the only change from System 1. The transmission 
was lowered at frequencies from 10 to 1000 eps. The second and 
third transmission peaks were split. 
Case 21. The System was the same as Case 20, except for the 
following: 

a. Mount stiffness decreased. 
The transmission over the freyuency range studied was lowered by 
approximately 10 db. 


Case 22. The System was the same as Case 20, except for the 
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following: 

ae Thiokol R. D. mount replaced Portsmouth mount. 
The first transmission peak was lowered; the second split peak was 
removed; and half of the third split peak was removed. 
Case 23. The System was the same as Case 20, except for the 
following: 

a. Mount damping increased. 

b. Mount stiffness decreased. 

ec. Hull damping increased twenty times. 
The first transmission peak was shifted below 10 cps, and 
transmission over the entire frequency range was lowered by 10 db. 
Case 24. The System was the same as Case 20, except for the 
following: 

ae Viscous foundation damping was added. 
The split transmission peak at 340 cps was lowered, amd the 
transmission from 1000 to 10,000 cps was lowered by 5 to 10 db. 
Case 25. The System was the same as Case 20, except for the 
following: 

a.» Thiokol R. D. mount replaced Portsmouth mount. 

b.-. Thiokol R. D. foundation damping was added. 
The first transmissdon peak was removed while the second and third 
peaks were lowered. The transmission was considerably reduced at 
frequencies above 300 cps. 
Case 26. The System was the same as Case 20, except for the 
following : 

ae Thiokol R. D. foundation damping was added. 

b. Mount stiffness decreased. 


The first transmission peak was shifted below 10 cps. With the 
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exception of the second peak, the vibration transmission was 
considerably reduced over the entire frequency range. 
Case 30. This was the basic case involving System 3. Foundation 
mass increase was the only change from System 2. The transmission 
was considerably lowered over the entire range except for peaks at 
800 and 4000 eps. 
Case 31. The System was the same as Case 30, except for the 
for Lowing t 

ae Mount stiffness decreased. 
The first transmission peak was removed, and the transmission 
reduced by 5 to 10 db. over the frequency range. 
Case 32. The System was the same as Case 30, except for the 
following: 

ae Thiokol R. D. mount replaced Portsmouth mount. 
The first, second, and third transmission peaks were removed. 
Case 33. The System was the same as Case 30, except for the 
Rellowilng:< 

a. Mount damping increased. 

b. Mount stiffness decreased. 

e. Hull damping increased twenty times. 
The results were the same as Case 31. 
Case .34. The System was the same as Case 30, except for the 
following: 

ae Viscous foundation damping was added. 
The first transmission peak was split; the peaks at 800 and 2500 
eps were lowered. The peak at 4000 cps was removed. 
Case 35. The System was the same as Case 30, except for the 


following: 
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ae Thiokol R. D. mount replaced Portsmouth mount. 

be. Thiokol R. D. foundation damping was added. 
With the exception of the transmission peak at 200 cps, the 
vibration transmission was considerably lowered over the range of 
frequencies. 
Case 36. The System was the same as Case 30, except for the 
following s 

ae Thiokol R. D. foundation damping was added. 

be Mount stiffness decreased. 
The first transmission peak was removed, and the transmission 
reduced by 5 to 10 db. up to 500 cps. Above 500 cps, the 
transmission was further reduced. 
Case 50. This was the basic case involving System 4. The 
machine and foundation masses were ten times the values in 
System 1. This System was used as the basis of comparison for 
larger machinery installations. 
Case 51. The System was the same as Case 50, except for the 
following: 

ae Mount stiffness decreased. 
The first transmission peak was shifted and lowered. The 
transmission was reduced by 5 to 10 db. over the 
Case 52. The System was the same as Case 50, except for the 
following: 

ae Thiokol R. D. mount replaced Portsmouth mount. 
The first transmission peak was removed, and the transmission was 
lowered by 5 to 10 db. over the entire range. 
Case 53. The System was the same as Case 50, except for the 


following: 
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ae Mount damping increased. 

b. Mount stiffness decreased. 

ec. Hull damping increased twenty times. 
The first transmission peak was considerably lowered. The 
transmission was lowered over the entire range by approximately 
10 db. 
Case 54. The System was the same as Case 50, except for the 
following: 

ae Viscous foundation damping was added. 
The first transmission peak was lowered; the third peak was 
eliminated. The transmission was lowered by 5 to 10 db. from 1000 
to 10,000 cps. 
Case 55. The System was the same as Case 50, except for the 
following? 

ae Thiokol R. D. mount replaced Portsmouth mount. 

be Thiokol R. D. foundation damping was added. 
The first transmission peak was removed, and the transmission was 
reduced by 10 to 15 db. to 300 cps. Above 300 cps, it was reduced 
by 20 to 40 db. 
Case 56. The System was the same as Case 50, except for the 
following: 

ae Thiokol R. D. foundation damping was added. 

be Mount stiffness decreased. 
The first transmission peak was lowered and shifted to 12 cps. 
The second peak was reduced by 35 db. At higher frequencies the 
transmission was extremely low. 
Case 57. The System was the same as Case 50, except for the 


following: 
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ae Thiokol R. D. foundation damping was added. 
be Mount stiffness decreased. 
ec. Mount damping increased. 
The first transmission peak was removed. The second peak was 
reduced by 35 db. At higher frequencies the transmission was 


extremely low. 
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CHAPTER IV 
DISCUSSION OF RESULTS 
4,1 Preliminary Discussion 

The most important single result of this thesis is the technique:’ 
of analytic description of a complex shipboard machinery installation. 
This technique formulates a workable analytic tool for the study of 
vibration transmission characteristics of these complex systems. In 
essence this technique is the four pole parameter method. It is recom- 
mended that this method of: analysis be considered in future design inves- 
tigations. 

4.2 Systematic Changes 
4.2.1 Mounting 

The properties of the resilient mounting had a profound effect 
on the vibration transmission. A study of the photo plots showed that 
stiffness in the Portsmouth mounting had a dominant iffect in the cases 
in which it was altered. When stiffness was increased, the transmission 
at all frequencies was increased; when stiffness was decreased, the 
transmission at all frequencies was lowered by 10 db. It is recommended 
that the design engineer make stiffness as low as possible and yet 
consistent with shock, alignment, and other factors influencing mount 
design. 

Damping in the Portsmouth mounting was not nearly as influential 
as mount stiffness. However, increased damping was effective in decreasing 
the transmission peaks at lower frequencies. It is recommended that an 
effort be made to increase: the damping of the Portsmouth mount by increas- 
ing the material damping or by introducing an external damping device in 
parallel with the mount. 


The replacement of the Portsmouth mounting by a synthetic rubber 
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mounting of Thiokol R. D. material showed interesting results. ‘The 
Thiokol R. D. mount lowered the first transmission peak and eliminated the 
second. Therefore, the use of the Thiokol R. D. mount appears desirable 
in systems with high noise transmission at low frequencies or in systems 
in which low frequency noise suppression is paramount. The stiffness of 
Thiokol R. D. increases exponentially with frequency (18). This decreases 
the»effectiveness of Thiokol R- D. in suppressing vibration transmission 
at higher frequencies. The development of a material to correct this 
undesirable trend of Thiokol R. D. is recommended. 
4.2.2 Foundation 

Of the various foundation parameters the mass most affected the 
vibration transmission. Three different foundation masses were investi- 
gated with the result that vibration transmission was lowered as the 
foundation mass increased. Within the limitations of weight and space, 
foundations of vibrating machinery should be made as massive as possible. 

Viscous foundation damping and damping utilizing the Thiokol R. D. 
material were studied. Foundation damping was effective at frequencies 
above 1000 eps for the viscous type and above 800 eps for the Thiokol R. D. 
type. However, the vibration transmission was reduced more by the Thiokol 
Re D. type. 

Foundation damping is recommended in installations where high 
frequency noise transmission is a problem. 
4.2.3 Hull 

The effects of increasing the hull damping first became apparent 
in the mathematical model used when the damping was increased by an order 
of twenty. As hull damping was further increased, the vibration transmission 
was reduced at all frequencies. 


These results are at best qualitative. In this model lumped 
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parameter damping was used, and no attempt was made to study the effects 
of distributed damping. In practice the latter would be the type applied 
to the hull. A study to quantitatively evaluate the effects of distributed 
hull damping would be useful in the extension of an analytic study. 
4.2.4 Combined Changes 

In a series of cases studied the Portsmouth mount stiffness was 
decreased, the mount damping increased, and hpll damping increased twenty 
times. This was done with the four systems studied. In general, as the 
masses in the systems were increased, the vibration transmission was 
reduced. The results of this series study were compared with those cases 
in which mount stiffness alone was changed. This comparison showed that 
the change: in mount stiffness overrode the changes in mount damping and 
hull damping. 

A single case of decreased Portsmouth mount stiffness and increased 
mount damping showed that the vibration transmission was lowered by exactly 
the same amount as the case in which stiffness alone was decreased. This 
further supported the conclusion that mount stiffness was one of the most 
important factors to be considered in any vibration transmission study. 

When a mount of Thiokol R. D. was placed in parallel with a 
Portsmouth mount, the transmission reduction properties were not improved. 
The addition of the Thiokol R. D. mount actually increased the transmission 
at frequencies greater than 1000 cps. This case was again studied with 
the Portsmouth mount stiffness decreased. Results showed that this case 
was not as good as that in which the single Portsmouth mount with decreased 
stiffness was examined. Therefore, it is concluded that placing the 
Thiokol R. D. mounting in parallel with a Portsmouth mounting would not be 
advantageous in reducing vibration transmission. If parallel mountings 


are to be used, as recommended by Snowden in Reference (18), care must be 
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used in the selection of the high damping synthetic rubber material. 

In another series of cases the Portsmouth mount was replaced by 
a Thiokol R. D. mount and Thiokol R. D. damping was applied to the founda- 
tion. This was done with the four systems studied. As the masses in the 
system were increased, the vibration transmissinn decreased markedly. In 
any one system the combination of Thiokol R. D. mount and Thiokol R. D. 
foundation damping was most effective in reducing transmission. It is 
recommended that this combination be investigated in an experimental 
installation. 

A Portsmouth mount with decreased stiffness and Thiokol R. D. 
foundation damping comprised the changes examined in another series of 
cases. Again this was done with the four systems. The vibration 
transmission was significantly decreased as the masses in the system were 
increased. Within a particular system this combination was less effective 
at frequencies of less than 40 cps than the Thiokol R. D. mount and Thiokol 
R. D. foundation damping combination. However, it was superior at 
frequencies higher than 40 eps. This combination should also be investi- 
gated in an experimental installation. 

A Portsmouth mount with decreased stiffness and increased damping 
and Thiokol R. D. foundation damping comprised the changes in the final 
series examined. The only transmission change:; from the previous series 
was the elimination of the first transmission peak. 

4.3 General Considerations 

The results of the systematic study largely verified existing 
concepts in reducing vibration transmission in shipboard machinery 
installations. This would indicate that the assumptions and the four pole 
parameter method of analysis were sufficient to establish the vibration 


trends and characteristics proposed at the beginning of this thesis. 
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As new ideas for reducing vibration transmission are introduced, it is 
proposed that they be investigated by means of the technique of this 
thesis. 

The four pole parameter technique enables use of empirical data 
obtained independently for each mechanical element in the vibration path. 
Empirical data can be used for elements which defy analytic description. 
It would appear that the use of empirical data for the element four poles 
would more closely approximate the vibration characteristics of an actual 
system. 

An interesting extension of this thesis would be an optimization 
study to determine the best combination of elements in the vibration 
transmission paths. Since shock specifications also govern the element 
selection, the optimization must include shpck as well as vibration 
reduction considerations. 

The construction of a model of the approximate system, as 
represented in Figure A=3 of Appendix A, would permit an experimental 


verification of the four pole parameter technique. 
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CHAPTER V 
CONCLUSIONS 
As logically developed in Chapter IV the following statements 
are listed as the conclusions of this thesis: 
1. The four pole parameter technique provides a sufficiently 
powerful method of analyzing vibration transmission. Furthermore, 
the four pole parameter technique is easily adapted to a digital 
computer program. 
2- For the purpose of vibration analysis a complex machinery 
installation can be represented by a relatively simplified system. 
This representation is: sufficient to establish vibration trends 
and characteristics. 
3. The stiffness of the resilient mounting has a dominant effect 
on the vibration transmission at all frequencies. 
4. The only effect of increased damping in the Portsmouth mount 
was to lower the first vibration transmission peak. 
5- A mount fabricated from a high damping synthetic rubber such 
as Thiokol R. D. is highly effective in reducing vibration 
transmission at low frequencies. The same mount is not as effective 
as the Portsmouth mount at higher frequencies. 
6. As the mass of the foundation is increased, the vibration 
transmission is Lowered. 
7+ Foundation damping was effective at frequencies above 1000 
eps for the viscous type and 800 eps for the Thiokol R. D. type. 
The Thiokol R. D. type was more effective. 
8. Increasing hull damping decreases vibration transmission at'all 
frequencies. 


9. Placing a Thiokol R. D. type mounting <in parallel with a 





CZ 
Portsmouth mount as suggested in Reference (18) actually had a 
detrimental effect on vibration reduction. 
10. Replacing a normal mount by a Thiokol R. D. mount and applying 
Thiokol R. D. damping to the foundation were most effective in 
reducing vibration transmission. 
li. Decreasing the stiffness of a Portsmouth mount and applying 
Thiokol R. D. damping 36 the foundation were superior to the 
system mentioned in the previous conelugion in reducing vibrations 
at frequencies above 0 eps. 
12. A Portsmouth mount with decreased stiffness and increased 
damping and Thiokol R. D. damping applied to the foundation were 


the most effective combination in overall transmission reduction. 
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CHAPTER VI 
RECOMMENDATIONS 
As stated or inferred in Chapter IV the following are the 
recommendations of this thesis: 
1. The four pole parameter method of analysis should be considered 
in future design investigations involving vibration transmission. 
2. As new ideas for reducing vibration transmission are introduced, 
they should be investigated by means of the technique of this thesis. 
3. Empirical data should be used for element four pole parameters 
to more closely approximate the vibration characteristics of an 
actual system. 
4, An optimization study to determine the best combination of 
elements in the vibration transmission paths should be conducted. 
It should include shock as well as vibration considerations. 
5. A model of the approximate system of this thesis should be 
constructed to permit experimental verification of the four pole 
parameter technique. 
6. Mount stiffness should be as low as possible and yet consis- 
tent with shock, alignment, and other factors influencing the 
design. 
7. The material damping of the Portsmouth mount should be 
increased, or an external damping device placed in parallel with 
the mount. 
8. The Thiokol R. D. mount should be used in systems with high 
noise transmission at low frequencies or in systems in which low 
frequency noise suppression is paramount. 
9. A high damping mounting material should be developed in 


which stiffness remains constant or decreases with frequency. 
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10. Within the limitations of weight and space, foundations of 
vibrating machinery should be as massive as possible. 
11. Foundation damping should be utilized in installations where 
high frequency noise transmission is a problem. 
12. A study to quantitatively evaluate the effects of distributed 
hull damping should be conducted. 
13. Caution should be exercised inthe selection of a high damping 
synthetic rubber mount to be placed in parallel with a standard 
mount 
14. The effects of replacing a normal mount by a Thiokol R. D. 
mount and applying Thiokol R. D. damping to the foundation should 
be investigated in an experimental installation. 
15. The effects of decreasing the stiffness of a Portsmouth mount 
and applying Thiokol R. D. damping to the foundation should be 
investigated in an experimental installation. 
16. The effects of decreasing the stiffness and increasing the 
damping of a Portsmouth mount and applying Thiokol R. D. damping 
to the foundation should be investigated in an experimental 


installation. 
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APPENDIX A 


DETATLS OF PROCEDURE 
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A-l The Beam-Mounting System 
A~l.1 Experiment 

This experiment was conducted to determine the isolation charac- 
teristics of a typical Admiralty mounting. The experimental system consisted 
of a fifty pound beam driven by a harmonic vibration generator. The beam 
was supported by the mounting which rested on a "rigid" foundation. The 
components and instrumentation are illustrated in Figure A-l. 

The foundation system was not, in fact, rigid. However, its natural 
frequency was much lower than the rigid body mode frequency of the beam~ : 
mounting system. This insured no resonance interaction, and the foundation 
MASS, I, in Figure A-1, responded in a rigid manner. The vibration 
generator imparted a velocity to the mass Mae The velocity V., into the 


3 


and the velocity V. at the free side of the beam 


mounting of compliance C 1 


ak 


were picked up by the barium titanate accelerometers, A, and AL respec- 


3 
tively. The isolation in this case is the ratio of V. to V.. Both this 


2 Ak 
and the driving frequency were recorded, and the isolation was converted 
to decibels. This was repeated over the frequency range of 10 to 3000 cps. 
To obtain a more realistic model in which the beam was driven 
from above, as a machine normally een! be, use of the transmissibility 
theorem was necessary. This stated that the transmission (or ESD ) 


in one direction for linear elastic systems is equivalent to the force 


transmission in the opposite direction: 


bp ed 
on 


F 
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This new system could be represented: 
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Figure A-1 
PRELIMINARY ExPERIMENT 
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The physical constants of the system of the experiment were as follows: 
Beam: Length: L = 60 in. 
dasekness: t = ./> in. 
Weight: Wa =) Ones 


Frequencies: 


w=, [EE 


where: E = 30 x 10° psi 


4 
I = moment of inertia. = — 
° 
@ = mass density 


Describing Functions: 


2 i) and A_ given in Reference (16) 
Ndi n 


= «Eo 
Mounting: Compliance: Cc. = I060 1b. 


A=-1.2 Theory 


The beam-mounting system was represented schematically: 
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[S, ] was the matrix of four pole parameters for the beam and 


L 
[S, ] for the mounting. 
Then: 
= 
vi “1 
at 


or by individual steps: 
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The matrices in. terms of their four poles were expressed: 


1. 
ea 
S = 
aE L., 
el 


1. 
i 

1. 
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and 


= 


: (4) 


As developed in Section 2.2.2, the general four pole parameter 


matrix is? 


1 b 1 fe 
=) ee 
om tp 
GQ = 
ral b 1 tT 
eel com x 


In this proof of the experiment the mounting was treated as a 


Simple spring. From Appendix A-2, the four poles for a simple spring are: 


s |= (5) 


v while the mobility is Jue, 


The beam was analyzed following the techniques of Ruzicka and 


The compliance of the spring is C 


Cavanaugh (4). With negligible compression assumed in the transverse 
direction and from the transmissibility theory, the four poles of the beam 


were evaluated :? 


iL ig 1 
(> =(2 =1 
F V 

b 
(M,, ) O 
yam ie 


oo = 


7 L 
ne M(>) 
where M(=) was the mobility at the center of the beam. This was taken 


from Reference (4): 
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The matrix [S, J for the beam was then written: 


i 
s, | = M3) (7) 
O 1 
Equations (5) and (7) were substituted in equation (1): 


1 





Fy Al uz) 1 O Te 
2 
= (8) 
jw) 
v4 O al j C, ul Vv. 
| =i eins 
To check the experimental results a plot of isolation I, which is 
i 
defined as S versus frequency was desired. Equation (8) was rewritten: 
3 
jwC V 
1 3 
F, = (1+ ) F =e 
aL L 3 L 
MZ) MS) 
Ge ( jw, ) F, - V3 
Station 3 in the experiment was a rigid support; therefore, V, = O. 


5 


The isolation was obtained: 





Bl JC 4 
id =, zit iL 
5 M(5) 


The mobility of the beam is expressed in Equation (6). 
A-1.3 Comparison of Experiment and Theory 

The curves of isolation in Figure A-2 show that the four pole 
parameter technique provided excellent theoretical verification of the 
experimentally measured values of isolation. These results indicated the 


value of the four pole analysis and suggested the use of it in the study 
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FIGURE A-2 
PRELIMINARY EXPERIMENT 
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of more complex systems. 
A-2 Development of Four Poles 


A-2.1 Mass Four Poles 





The performance equations are 


V, = Vp Gy 


Fo- Fo =m een (2) 


from the application of Newton's Law. 
pinee the applied foree is usually considered as harmonic, the 


Other force and the velocities will also be harmonic functions. 


a jut m jwt 
FL = Fio e F, = Fo e 
e jut Ss jut 
Lemos ee) 
Figo» E50? Vio» and V.. are time independent but may be complex. 


20 


Substituting the expression for V,, Equation (2) becomes: 


kee 


F F_ + (jwm) V 


i © 2 


while (1) remains 


a = © 
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By inspection the four pole parameters are: 


€e., =1l, € 


11 1 


A-2.2 Spring Four Poles 


mail 


The performance equations are 





FL = F, (1) 


because the spring possesses no inertia and 


Fy = k [x - Xp where x is displacement. (2) 


The velocity is the time derivative of the displacement. 


From the harmonic velocity expressions of the previous section: 


7 just = jwt 
Vi = Vio e V5 = V50 2 
and 
Z jut = jut 
x, = (vy, e at x, -(v,, e ag 
or 
oe 
1 jw 2° jw 


Equations (1) and (2) are rewritten: 


P= 
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A=-2 .3 Damper Four Poles 





The performance equations are 


PL =F, 


because the damper possesses no inertia and 


from the definition of damping. 
These are rewritten 


a 


va = 


A-3 Parallel Matrix Combination 


If n systems are connected in parallel then a single four pole 


(1) 
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matrix can be written connecting the common input and output points. If 


this matrix is written Is, then 


— 


A AC 

5 aaa 
S = (1) 
° a 

B B 


Let the r“” matrix, [s], be given by: 
S| = (2) 


where r varies from 1 ton. Then, in (1) above 


a=) | (3) 








B=) 1 (14) 
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A-4.1 The Machine 

As stated in Section 2.4.1 the machine was represented by a 
uniform elastic beam driven by a harmonic force generator. This system 
was selected because it was multi-resonant and possessed distributed 
parameters. 

The four poles of this system are those briefly described in 
Appendix A-1.2 for the beam in the preliminary experiment. These four 
= a = QO, and eé 


oe aaa 
M(S) 


poles were E44 = 1, E10 p2 = L. From the method 
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of Reference (4), the beam mobility M3) will be derived. 
The basic equation of free vibration is: 
4. 92y 
Hate 5 = 0 (1) 
Ox ot 
where: I = moment of inertia of beam's cross section 
K = modulus of elasticity 


4 = mass density 


> 
Il 


beam cross-section area 
y = transverse displacement 
x = distance along length of beam 


L = length of beam 


I 


i time 


Assuming a harmonic solution: 


y = t(x)e”* (2) 


If this is substituted, equation (1) becomes: 


where 


2 
4 Aw 
Be =P 
and W= frequency of vibration. 


The general solution of (3) is: 


f(x) = C, (sinh @x + sin 4x) + C, (cosha x + cosgx) + C., (sinh 


1 3 


Bx > sin 6x) + C), (cosh 2 x — coss x) (4) 


The boundary conditions for a free-free beam are: 


ae O 


ax 


do (5) 
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(8) 


| 
© 


a 


Equations (5) and (6) require no shear forces at the ends of the beam; 
equations (7) and (8) state that there are no end moments. 

When the boundary conditions are applied to equation (4), the 
characteristic frequency equation is obtained. 


cos B. L cosh ea liye) a! (9) 


This is a transcendental equation. Corresponding to each eigenvalue B. 
will be a value of f(x) . The term f(x) represents the mode shape in 
which the beam vibrates at the frequency defined by Pr’ The relationship 


between £ Ax) and 8, is: 


ffx) = C,[cosh B, x + cos 8, — x ( sinh Z, x + sing, x) | (10) 


where 
sinh £. L + sin 2, L 
ol 7 cosh £ L- cos 6B. L ee) 
C, = arbitrary constant 
Define :-: , 
f(x) = - d(x) ° 


Then in equation (10): 
(x) = cosh B x + cosB x ~& (sinh 2) x + sin 8 x) 


The term @ (x) will satisfy the equation: 
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y 
ao. 


1 in Bn On 
It is a normal function which forms an orthogonal set in the interval 


©sxtL. The condition of orthogonality: 


S44, dia —20 m ¢ n 


However, when n = m;: 


L 
Son) & =F 1G, DP 


D. Young and R. Felgar have evaluated these ® functions for beams with a 


| 


variety of boundary conditions (16). 

With the values of O. and the corresponding values of WwW the 
forced vibration motion is considered. 

If a harmonic force Fe" is applied to the beam, the displacement 


according to reference (1) is: 


an 
y(x) = ry bn) rigid body motion C2) 
= lw* Swe 
Jwt 


The rigid body motion is f) (x) e The velocity is the derivative of 


the displacement: 


: 2) 

2 jwFk aulce - rx . 

ee a > BB t Why (x) eM (13) 
n= We 
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Prem thew@efinition of mobility at the driving point p: 





F(p jut 
Foe 
“(p) 
» bot 
yy re 
B nelW, -w O 


As a rigid body the beam acts like a mass and possesses the mobility of 


this lumped parameter. From Appendix A-2.1: 


“J 
mw 


Therefore, in equation (14): 


Mass Mobility = 


Jui, _ 3 
By mW 





The final expression for beam mobility is: 


= 2 
; ) ; 
Mp) = 24 a Gui?) ng (15) 


Since the beam in the machine representation was considered to be 


dimiven at its center, the four pole ¢« required the inverse mobility 


12 


at the center of the beam. This was expressed from equation (15) as: 


L, gw ) o2,(5) j 
MS) = an 22 mw 
"B n = LW —w m3 

The machine could now be represented in terms of its four poles. 
A-4.2 The Mounting 

According to Section 2.4.2, two different mountings were studied. 
The first was the standard Portsmouth mounting which was analyzed following 
the method of Wright and Hagg (17). 


The empirical results of tests on actual mountings were used by 


Wright and Hagg. The mount characteristics were approximated by use of the 
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electrical impedance analogy, one used extensively in vibration analysis 
(8). The mounting was represented by a shunt impedance consisting of real 
and imaginary parts. In the language of the analogy, impedance is defined 


as a force to velocity ratio. 


MN 
I 
<j] hy 


(1) 


In this case input impedance was considered for the approximation. 
From the analogy mechanical damping corresponds to electrical resistance 
or the real part of the impedance. Mass and spring-like properties 
correspond to electrical reactance or the imaginary part of the impedance. 
Thus the approximation included all the important lumped parameters. 

These parameters, damping constant R and spring stiffness K, were 


established to be: 


K 
R= ontQ (2) 
0 
mre (a 
0 fo 
K= re (3) 
0 
g(1 + 2) 
t 
where ; 
fy = rated natural frequency of the mount at its rated Load, cps, 


empirically determined. 


Q. = maximum transmissibility at f., dimensionless, empirically 


0? 
determined. 


f, = frequency of the first high frequency transmission peak, cps, 


t 
empirically determined. 
-* rated mass load on the mount, lbm, empirically determined. 
f -= driving or exciting frequency, cps. 


The shunt impedance Z is therefore given by: 





ee) 


(4) 


The impedance analogy indicates the following equations hold for 


a shunt impedance? 


eat 
ee, ? 





ian - eae 





Open Cinevit 


pes Oine @CacenaG 


Equations (5) to (7) are translated into four pole parameter 


notatbhon as in Section 2.2: 


a we oe 

a Bae 2 
(M,.) 
aes 


station 2 free 


(5) 


(6) 


re 


(8) 


Sy) 


(10) 


Therefore, the four pole parameters for this standard Portsmouth 


mounting are 


Combining equations (2), (3), and (4) the following expression for e¢ 


is obtained. 


vag 
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Way lL + JQ, 


i ae a) 


a 

Q +1 

where wis the forcing frequency, rad/sec. 
The second mounting considered was treated similarly except that 


different expréssions were used for R and K. These were derived by A. O. 


Sykes (14): 


2rG S 
R = S20 ____ (3) (12) 
wW({l-r) 
S 
K = (2) (13) 
WNere ? 
5 = mounting cross-sectional area 


lL = mounting length 


G2 
It 


mounting elastic modulus 


K 
It 


mounting damping factor 

w= forcing frequency 

The elastic modulus and damping factor were given as empirical 
results for a synthetic high damping rubber called Thiokol R. D. by J. C. 
Snowden (18). 

Approximate formulae were determined for G and r from the empirical 


datas 
G = 72 .6ur** (14) 


= Cem (15) 
Thus R and K could be represented as functions of frequency and the new four 


pole parameter E54 evaluated. 


a ae 
“~  a(2) (ae + 1) 





(54) 
where 


or 


Ll - a 


Ac 


A-4.3 The Foundation 

As stated in Section 2.4.3, the foundation was considered a beam 
resting on two columnar structures. The two columns presented two possible, 
identical paths for vibration transmission. ‘Since the paths were identical, 
only one was studied explicitly. 

The four poles of the foundation structure were a series combination 


of the individual four poles of a beam and column. 


dl: 
al S75 al M(L) al O 
Coy Eno 0 ab Ma aL 


The combined element four poles are 


M 

€4, = 1 EP E15 = Ty fo. = Ma and Coo = Ls 

The term M(L) is the transfer mobility of the beam relating the 
driving point with the junction of the beam and column structure. Initially, 
this beam was considered as end-supported for the purposes of deriving its 
mobility. Correspondence with the Admiralty Engineering Laboratory, where 
an experiment driving an end-supported beam was being conducted, indicated 
that for vibration analysis he end-supported beam behaved more like a 
free-free beam. For this reason the foundation beam was treated much as 
the beam in Appendix A-+.1. The expression for beam mobility is: 


M(L) “~~ Zz A Nae 


ne lw -w Hoel 


where the d's are the normal functions at the driven point (2) and the beam 


column junction (L). Again these functions were taken from Reference (16). 
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From Reference (1) the expression for the lateral deflection of 


a column being driven by a harmonic force Reo 1s 


Pere : 
P = mass density 
L= column length 
A= colwmn crosss;section area 
The column mobility can be derived using the same methods used for the 


beam in Appendix A-4.1. The velocity is: 


igi ib 

‘ JU Nee ——— 

=) ata <1 a(-1)° 
y(x) = Jwfe + = 5 5 
fe n=1 (We -W~ ) 


Now the definitions of mobility and the rigid mass mobility are applied. 
From this the mobility of the colum is obtained: 
zs n=1 
; : a 
My =i * pia )_ ee 
e ) 
A~4.4 The Hull 

As stated in Section 2.4.4, the work of D. V. Wright was consulted 
again for a workable analytic representation of the ship's hull (7). 

Wright considered the hull a series impedance, which he represented 
as a result of analyzing empirical data obtained from typical hull impedance, 
measurements (11, 20). As in the case of the mounting the impedance 
consisted of a real and imaginary part, the resistive and reactive properties. 
These were represented as complicated quantities dependent on nine separate 


lumped parameters. 


PS 
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where Ro» GR eesmnN wom sn Ke iy Mh are the nine lumped parameters and 


2 Spams. 9D 


f is the driving frequency. 


1? 


The nine parameters were considered circuit elements in the 


impedance analogy, making up the series impedance Z. 





The impedance analogy gives: 


Fin 
ht (4) 
in 
Van ~ Vout (5) 
Fe = Foyp? Open circuit (6) 


Equations (4) through (6) are translated into four pole parameter 
notation by the method of Section 2.23 


Vv, = V, (7) 


1 E, station 2 blocked (8) 


7) 
i 





A~4.5 


ai 





(My) Ne 


The 


The 


(5%) 


F 


-73-2 (9) 


four pole parameters therefore become: 


=1 (Ly) 
=% (11) 
aye (12) 
27 (13) 
Water 


As mentioned in Section 2.4.5 the water, considered only as a 


resistance to the propagation of acoustic waves, has a specific resistance 


R* dependent on the water density and propagation velocity. 


This can be 


written, (2): 


where 


R¥ = 


R* 


we) 
I 


I> 
H 


(A 


Cc 


il 


| 
b> | bd 
1 
Y 
QQ 


(1) 


specific resistance 
Resistance 
Effective radiating area 
density 


velocity of sound in water 


Treating the water as a dashpot-type resistance (velocity 


proportional), the four pole parameters were obtained in Appendix A-2.3. 


(2) 


= 1 


(3) 


=) 


I 
bd] 


a (5) 
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R was evaluated by two methods. In one an effective area A was 
assumed cOnstant with frequency. 
R qa (6) 
The hull was considered also ag. a vibrating dtaphragm. The radius 


r of the effective vibrating area is given by: 
Cc 
ra (7) 


where w is the radiating frequency. 


Therefore: 


ae (8) 


Substituting (6) or (8) into (4) gave the desired four pole para~ 
MELEYrs. 
A-5.1 General Description 

Four variations of the general system of Figure A~3 were analyzed 
in this study. These variations were made in the relative sizes of the 
machine and foundation structures to minimize any scale effect in the 
trends of vibration transmission. Within each of these systems, element 
properties were varied. In addition, isolation elements were inserted at 
various locations in the systems. From these variations and insertions, 
trends to provide a guide for preferred design practice in the reduction 
of vibration transmission were established. 
A-5.2 Syatem 1 

With reference to Figure A-3 the properties of System 1 were: 

1. Machine 


Weight of beam, We = 50. Oo ia 


Mass of beam, m, = .1295 lbm 


B 


ee 





(59) 


FIGURE A-5 
GENERAL SYSTEM 
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Length of beam, L = 60.0 in 
Thickness of beam, t = ./> in 
Mounting 


Rated natural frequency of the mount at rated load, 


Maximum transmission at fy Qo = 20 


Frequency of the first high frequency transmission peak, 


: te ages 
eT a 


Rated load of mount, W = 50.0 lbf 
Foundation 


Weight of beam, W., = 30.0 lbf 
Mass of beam, m, = .O78 lom 


Length of beam, L = 36.0 in 
Thickness of beam, t = .75 in 
Column 

Cirewular Cross—-sece1con 

Length of column, L = 6.0 in 


Mass of column, m= -039 lLbm 
Weight of colum, we = Cas 


Huda 


ihe nine parameters are; 


R_ = 50 lbf-sec/in 
R, = -O1 lbf-sec/in 
R._ = 5000 lbf-~sec/in 


meray yi 


i | 





A-5.3 


A-5.4 


aD 
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X, = 500 lbf-sec/in 
W. = 10.0 lof 

Gane. /. 
K, = 10 16f/in 
W. = 15 lbf 
Pp 

6 
K, —eeseyO) 1bt/in 
6. Water 


R = 37,600 lbf~sec/in (constant with frequency) 
system 2 
With reference to Figure A-3 the properties of System 2 were: 
identical to System 1 excepe: 

Foundation 


100 lbf 


x 
tt 


259 lbm 


= 
I 


SOO alia 


Il 


t 


I 


i Sedan 

oystem 3 

With reference to Figure A-3 the properties of System 3 were: 
idenclegietersystem 1 except: 


Foundation 


W 


2 200miiba 


P5410 ibm 


ue 


L = 60.0 8in 


t = 3.00 in 
System 4 


With reference to Figure A-3 the properties of System 4 were: 
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Identical to System 1 except: 


Machine 

Wa =—epoo br 
m, = 1.295 lbm 
L = 60.0 in 

1 53 omen 
Foundation 

We = 300 bf 
mi, = .78 lbm 
L = 60.0 in 
ee: salle 


A-5.6 Element Insertions 
A-5.6.1 Thiokol R. D. Mounting 
In the systematic study the Thiokol R. D. mounting was either 
substituted for the Portsmouth mounting or placed in parallel with it. The 
properties of the Thiokol R. D. mounting were: 
Mounting cross sectional area, S = 9.43 oe 
Mounting length, L = 4.0 in 
A-5.6.2 Foundation Damping 
Viscous foundation damping which was constant with frequency was 


first considered. The damper was inserted between the foundation and the 


colum. The property of the damper was: 


System .1: 
Damping constant, R = 390 aeeasee 
System 2: 
lbf-see 
wigio7* in 
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system 3: 
a = 1560 1weasee 
sig 
System 4; 
2c ig Sees 
in 


Damping at the same position with Thiokol R. D. as the damping 
material was next considered. The dimensions of the Thiokol R. D. were: 


Cross sectional area, S = 20.0 ine 


iimemess, i°— 1.0 dn 
The thickness of the damping material was such that its compliance 
property could be neglected. 
A-5.7 pystematic Changes 
The systematic changes in the structural foundation path are best 
followed by an examination of the individual cases. 
Case 1. System 1 used. This was the basic case for System 1 with 
which the cases with property changes could be compared. 
Case 2. System 1 used. Mount damping was increased by decreasing 
Ro Ome sr 
Case 3. System 1 used. Mount stiffness was increased by setting 
W. = 201.062 and Wi = 12,566.384. 
Case 4. System 1 used. Mount stiffness was decreased by setting 
WJ. = 50.266 and Wy = Byers 596. 
Case 5. System 1 used. ‘The Thiokol R. D. mounting replaced the 
Portsmouth mounting. 
Case 6. System 1 used. The hull damping was increased by setting 
R. = NG@Ot7 nh = 02, ander. — 2000. 


a 2 


Case (. System 1 used. The mount damping was increased by setting 
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and W, = 3,141.596. This case was thus a combination of cases 2 and 


= 2.5, and the mount stiffness was decreased by settin WW, = 50.266 


is 
Case 8. System 1 used. The Thiokol R. D. mounting was used, and 


the hull damping was increased by setting R. = LOO hear 7 O2eand 


af 


R, = 1000. This case was thus a combination of cases 5 and 6. 
Case 9. System 1 used. Mount damping was increased, mount stiffness 
decreased, and hull damping increased. This was accomplished by 
setting Q, = 2.5, Ld, = 50.266, W = 3.141.596, R, = 1000, R, = 1.0, 
and R, =w20)000; 
Case 10. System 1 used. Viscous foundation damping was inserted 
in the system with R = 390. 
Case 11. System 1 used. The Thiokol R. D. and the Portsmouth 
mounting were placed in parallel. 
Case le. System 1 used. The dissipative resistance of the water 
was varied with frequency. See equation (8) in Appendix A-4.5. 
Case 13. System 1 used. The mount stiffness was decreased by i 
setting ud = 50.266 and L4 = 3,141.596, and viscous foundation 
damping was inserted in the system with R = 390. This case was 
thus a combination of cases 4 and 10. 
Case 14. System 1 used. Foundation damping with Thiokol R. D. 
was inserted in the system. 
Case 15. System 1 used. The Thiokol R. D. and the Portsmouth 
mounting were placed in parallel. The stiffness of the Portsmouth 
mounting was decreased by setting we = 50.266 and W, = 3,141.596. 
This case was thus a combination of cases 4 and 11. 
Case 16. System 1 used. Hull damping was increased by setting 

6 


= 100, and R, =2x10. 


R_ = 100,000, R A 


AL 
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Case 17. System 1 used. Foundation damping with Thiokol R. D. 
was inserted, and the Thiokol R. D. mounting replaced the 
Portsmouth mounting. This case was thus a combination of cases 
5 and 14. 
Case 18. System 1 used. Foundation damping with Thiokol R. D. 
was inserted, and mount stiffness was decreased by setting 
uy) = 50.266 and ud = 3,141.596. This case was thus a combination 
of cases 4 and 1}. 
Case 19. System 1 used. Foundation damping with Thiokol R. D. 
was inserted, mount stiffness was decreased by setting us = 50 .266 
and WJ, = 3,141.596, and mount damping was increased by setting 


Q 


ae 2.5. This case was thus a combination of cases 19 and 2. 


Case 20. System 2 used. This was the basic case for System 2 
with which the cases with property changes could be compared. 

Case Gl. System 2 used. Mount stiffness was decreased by setting 
ud = 50.266 and ly = 3,141.596. 

Case 22. System ec used. The Thiokol R. D. mounting replaced 

the Portsmouth mounting. 

Case 23. Systeme used. Mount damping was increased, mount 
stiffness decreased, and hull damping increased. This was 
accomplished by setting Q = 2.5, up = 50.266, Wy = 3,141.596, 

= i. .0, ana kn 


e 


Case 24. System 2 used. Viscous foundation damping was inserted 


R_ = 1000, R = 20,000. 


i 
in the system with R = 390. 

Case 25. System 2 used- Foundation damping with Thiokol R. D. 

was inserted, and the Thiokol R. D. mounting replaced the Portsmouth 
mounting. 


Case 26. System 2 used. Foundation damping with Thiokol R. D. 
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was inserted, and mount stiffness was decreased by setting 
Us = 50.266 and Wy = 3, 141.596. 
Case 30. System 3 used. This was the basic case for System 3 
with which the cases with property changes could be compared. 
Case 31. System 3 used. Mount stiffness was decreased by setting 
Wy = 50.266 and W, = 3, UHI SG. 
Case 32. System 3 used. The Thiokol R. D. mounting replaced the 
Portsmouth mounting. 
Case 33. System 3 used. Mount damping was increased, mount 
stiffness decreased, and hull damping increased. This was 
accomplished by setting Q, = 2.5, Wh = 50.266, u) = Spee. 596. 


R, = 1000, R, = 1.0, and R, = 20,000. 


aL 
Case 34. System 3 used. Viscous foundation damping was inserted 

in the system with R = 1560. 

Case 35. System 3 used. Foundation damping with Thiokol R. D. 

was inserted, and the Thiokol R. D. mounting replaced the Portsmouth 
mounting. 

Case 36. System 3 used. Foundation damping with Thiokol R. D. 

was inserted, and mount stiffness was decreased by setting 

uy = 50 .266 and U4, = 3,141.596. 

Case 50. System 4 used. This was the basic case for System 4 

with which the cases with property changes could be compardd. 

Case 51. System 4 used. Mount stiffness was decreased by setting 

Ud = 50.266 and tW = 3,141.596. 

Case 52. System 4 used. The Thiokol R. D. mounting replaced the 
Portsmouth mounting. 

Case 53. System 4 used. Mount damping was @ncreased, mount stiffness 


decreased, and hull damping increased. ‘This was accomplished by 





(67) 


setting % = 265) UW = 50.266, Wi Bp lle 590; R. = 1000, R. = 1.0, 


at 
and. R, = 20.000). 

Case 54. System 4 used. Viscous foundation damping was inserted 

in the system with R = 1760. 

Case 55. System 4 used. Foundation damping with Thiokol R. i. 

was inserted, and the Thiokol R. D. mounting replaced the Portsmouth 
mounting. 

Case 56. System 4 used. Foundation damping with Thiokol R. D. 

was inserted, and mount stiffness was decreased by setting 


up = 50.266 and w = 3,141.596. 


O iv 
Case 56. System 4 used. Foundation damping with Thiokol R. D. 


was inserted, mount stiffness was decreased by setting W = 50.266 


and WW = 3,141.596, and mount damping was increased by setting 


Q. = 2.5. 


0 


A~6 Sample Calculations 


As noted in Section 2.5.2, the computational effort involved in 


the thesis was so great as to render it impossible without the use of a 


digital computer. However, the computational effort was largely needed 


in iterative calculations. To illustrate the complexity and length of 


the calculations involved in solving this problem, a sample calculation 


for one value of transmisston is presented here. 


For example, Case 1 properties are used with) the driving fvequenc,, 


at 10.3 cps or 65 rad/sec. 


The unsimplified four pole matrix equation for Case 1 follows: 





é dh. 

V: 0 al — 21 M a 

al Z, c 

3 J. 
Machine Mounting Foundation 
ili Ley ai O iY 
(1) 
1 

O 1 R iL Ue 


nb k aL ! Water 
ime toliowing are the properties for this case: 


1. Machine 


We = 10) sll loni 
fee= 60.0 in 
(Geel Oe celts al 


e. Mounting 


fy = 16 eps 
Q = 7.0 

i = lLOOO cps 
Wo 5OrO 


3. Foundation 


We = 30.0 lbf 
mn, = O78 lbm 
E = 36.0 10 


tee 2.75 in 
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uh. Colwm 


lL = 6-0 an 

ma =e 40 Lom 

W. =) AED Olen 
Dew Hult 


The nine parameters are 


R. = 50 lbf-sec/in 
R, = -O1 lbf-sec/in 
Ry = 5000 1lbf-sec/in 
xX. = 500 lbf-sec/in 
W, = 10.0 lof 

6 ; 
K. = ko lof/in 
W. = 15 lvof 
Pp 

6 , 
K, = 2 x 10 lpf/in 
6. Water 


R = 37,600 lbf-sec/in 
These properties are also listed in Appendix A-5.2. 

The individual four pole parameters must be evaluated before 
combination of the matrices. 
Machine: 

From Appendix A~4.1: 


oy ee 5 bate) aes 


M( 
Sh a a wh Su ee 


| 


where: 





@-a/ = 
n n mM, 


L 
QD (2) = COs = 


For W= 65 


L 
2 a 


sinh 4 L + sin ae 


cosh BL - cos 6 1 


and ob (2) are listed in Reference (16). 


I] 


1.21565 


= 751 a = 5080 o, =a 
= 1475 ae = 6780 o, = +1.42238 
= 2hho ud = #738 >, = 0 


ll 


= 3630 w= 19999) 5 = ~1-41366 


» using the above and listed properties 
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L ; 
M(S) =—- .1074 3 
il. ; 
lt, 7 ¢ 
M(3) 
Mounting : 


ay 


M 


om Appendix A-4.2: 


K 


WQ 1 + jQ 
a: O 
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Oo 
Qo + 1 


a ilar 


oe is calculated. 
P= 0 
Se ahs, 


Py = 0 


0, Se iihiae 


M9 = 0 


aSF 
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_@ 
1295 2, Ww _ 205 (65) (eHSoE) 
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For W= 65: 


2 105 
- (alee) 100 Ge) [ies 583.2 
7 ales n( 100.53) 
6283.2 


K # 1285 
i CO er lL + 73 
he ~ 1285 t9 + 1! 
_ = .00708 + .0h96j 
M 

Foundation: 


From Appendix A-4.3:; 


oO aly 
nny) = 0 ) QOH gy 
eB a= 1lWo -w mB 





The values of We and Q (1) were calculated as for the machine. The values 


of (2) are the same as for the machine. 


We 758 a = 10100 d, = +2.0 > = 92.0 
UZ, = 2085 ln = 14080 0, = -2.0 0. = +2.0 
W = O90 bd = 18780 ), = +2.0 Ys = -2.0 
Wy, = 6770 Wh = 24200 Q, = -2.0 Wy = +2.0 
W, = 10100 W,. = 30100 Qj, = 42.0 W,, = 32-0 
For W= 65: 
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146.704 


Od 
I 


_ dw _ J 
ae arty F) m uJ 
From Reference (1): 
ae 

ane an ie F 

W = 3.33 x 10? 

1 

os = 9.99 x 10° 

W, = 16.65 x 10° 
UW) = 23.31 x 10° 

a = 29.97 x 10? 

For W = 65 

2 n-1 
, - 365) > aCe 
= Ey 1 (we - 4095) 202 

We a7 0 

From Appendix A-4.4 
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A) Ep siele 65 
X = 500 + (5 0) le Pee UG, 
386 5 a 2 386 65 
+ (500) (eS NIBY or a8 
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X = ~ 14,898.3 
Z = 146.7 -— 14,898.3 5 


Water: 
From Appendix A=4.5 


R = 37,600 


Substituting the computed four pole parameters in equation (1) 


in polar form: 


| 13 I 
Fy il 9.30/35 ih O 2.96/ 0 NnOeT = 
ee 
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This equation will reduce to the form: 


(3) 


HH 


< 
a 
am 
a 


With reference to Figure V: 


v= 0 (4) 


From the four poles of the element representing the water: 
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In section 2.5.1 the transmission T was defined as? 


*5 
T = 20 108.6 - 
ab 


Therefore: 


Ale 


ill 

20 log —— 
Ome 

WAL 


From obtaining equation (3): 


Sa, Oe 
At Wes 65: 
T= + 14.85 


A-7.-1 Computer Logic 

The logical procedure comprising the program written to solve this 
problem can best be illustrated by the computer flow diagram in Figure A-4. 
A-7.2 scope 

The cathode ray tube oscilloscope is a convenient computer feature 
which produces a graphical representation of the problem answers. ‘The 
computer routine used to program this feature is not included in the 
program listing. Writtenhnin FAP, this routine projected a three cycle 
semi-Log grid on the scope, labeled the grid, and plotted the answers. 
Photographs, 35 mm. in size, were taken of the scope presentations. These 
photographs appear as results in Appendix B. 


A~7.3 Description of Subprograms 
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Figure A-4 


value of 
a 










Compute 
Machine 
4 Poles 





Compute 
Mounting 
4 Poles 


ete ee 


Compute 
Foundation 
h Poles 


Compute 
Foundation 
4 Poles 


Compute 
Water 
h Poles 





Compute 
Answers 


Increment 
OG) 


Computer Logical Flow Diagram 


is 





value? 







Endwot viob 





(76) 
1. MAIN This serves as the control program calling all subprograms 
when needed in the solution. 
2. RDIN This subprogram reads in the data including the system 
properties. It also calls the routine to place the semi-log grid 
on the scope. 
3. MAT! 1 This program computes the four poles for the first element 
in the system. In this case the four pole parameters for the 
machine representation were obtained. 
4, MAT 2 This serves the same purpose as MAT 1 for the second 
element or mounting in this case. 
>- MAT 3 This serves the same purpose as MAT 1 for the third 
element or foundation in this case. 
6. MAT 4 This serves the same purpose as MAT 1 for the fourth 
element or hull in this case. 
{- MAT 5 This serves the same purpose as MAT 1 for the fifth 
element or water in this case. | 
8. TRANS This program takes all the four pole parameters computed 
by the various MAT routines and stores them for use in the matrix 
multiplication. It sets up the matrices prior to multiplication. 
9. MAMUL This program reduces the complex matrices to one 
resultant matrix. 
1O. ANSW This program selects the matrix elements to obtain the 
answer from the final matrix of MAMUL. It also provides the routine 
for plotting the points on the scope. 
ll. CHOMEG This program tests the value of frequency and trans- 
mission to insure that the scope points are numerous enough to give 
an adequate plot. It also changes w before returning to the start 


of the program. 
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12. PRINO This program prints the values of frequency and 
transmission to conclude the program. 

A-7.4 Program Listing 
The computer program listing is included to facilitate its possible 


future use in design investigations. 





210 
240 


st ie) 
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MAIN PRUOGKAM 

DIMENSION Al25) 956425) 9C (20) 201 2U) 8b (400) 9G(40U) 9G1(40U) »V(20U) 9Vb 
1(20) sAA(2ZUU) sthis(2UU) »CC(2UU) sUDI1UU) sEE (LOU) 

CALL RUIN (ApBeCeVsAAsnbesCLsbiett) 

Gil) = AY) 

DO 210 L#l»s400 

MATRIX CALL SEQUENCE 


CALL MATL (GsAsBeCeVoLl sAAsbbesCCeDVett) 
CALL MAT2 (GeAsbBeCeVoLl sAAsbieCCeUDett ) 
CALL MAT3 (GeAsBeCeVoel sAAstbsCCsLUeEkt) 
CALL MATS (GsAsbeleVoel sAAsbieCCeUUeEkt) 
CALL MAID (GeAsbseCeVol sAAsbiaCC eUiUstt) 
TRANSLATION SEQUENCE 

CALL TRANS (VetE) 

MATRIX MULTIPLICATION 

CALL MAMUL (1eJeKoVeA) 

ANSWERS 

CALL ANSW (AsbLleGot eLoV) 

CALL FREQUENCY CHANGE 

CALL CHOMEG (Eol»sGeA) 

IF (G(L)7A(8)) 21092402240 

CONTINUE 

CALL PRINU (GlsGeteL) 

CALL EXIT 

END 

SUBROUTINE TELE 

CALL FRAME 

CALL GRID 

CALL GRID 

CALL GRID 

CALL GRID 

CALL GRID 

CALL GRID 

CALL TiTLe 

CALL VERT 

CALL HOR 

RETURN 

END 


SUBROUTINE KDIN (AseCesVUsAAsbisCCsDiecEt ) 
DIMENSION Al2d) ebl 2d) eC (20) sD ZU) sAAI CUU) shits (ZOU) »CC(20U) sUDI LULU) » 
lEE( 100) 
READ 213 e((ACT) oL=lL9lU) oS ( BEL) sl=le20) 9(C(L) st=1920) 9(D( 1) sl=lecV) » 
L(DO(1)sT=1210)) 
FORMAT (>SE13e/) 
CALL TELE 
RETURN 

END 


SUBROUTINE MATL (GsesAsbsCeUoel sAAsvniieLCaLvU sEL) 

BEAM MATRIX 

DIMENSION GO14U0U0) 8A 2d) sul ed) eClceVU) sUI CU) PAAI ZUYU) ebb (200) eCC(guUu) »s 
LDVU( 100) see (100) sVB( 20) 
V43=zleO0/(b(2)#G(L.)) 
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208 


209 
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DO 205 I[=lse 10 

IF(DCT) **2-G(L)**2) 20192022201 
VB(TPIRG(LIFC(LI##2/ (001) #¥2-U(L) Ft?) 
GO TY 205 

IF (G(L)#C C1) #*¥#2-Q500) 20492039204 
VB(1)=0¢4«0 

GO TU 205 

VB(IT)FA( Ss) 

CONTINUE 
V7BzVE(1L)tVdb(2)4+V¥b(3)4+V6(46)+V09)4+Vb(60)4+V0(07)+Vb(56) +V6(9)4+VB(10) 
V/I9=VIB/B( 2) 

V229=VI9—-V43 

TF(V229-OeV0) cUb6s208920/ 
V329=A(6) 

GO TU 209 

V329=A(2) 

GO TU 209 

V329=000 

V328=A(7) 

GO TO 212 

V328=1leO/ABSFIVZ229) 

EE( 1) V32u 

EEC } V329 

RETURN 

END 


not 


THIS IS MOUNT 
SUBROUTINE MAT2 (GseAsBeCeDoel sAAsBBeoCCesDDoEE ) 
DIMENSION 61400) 9A(25) 9KB(25) 9C(20) 9D( 20) 2AA( 200) 9BB( 200) »CC(200) » 


100(100) sEE (100) 


VA = (B14) **2%Kb(5))%*((Bl6)4+6(L) *A(1))7(B16)4B8(14)*A(1))) 
RE = (G(L)*8(3))/(VA*(B(3) #¥2+1¢0)) 

VE = RE*B(3) 

EE(2) = SOQRTF(VE**2+RE**2 ) 

EE(4) = ATANF(VE/RE) 

RETUKN 

END 


SUBROUTINE MAT2 (GeAsn eC evel sAAsdts es CCsDUDoE“L) 
DIMENSION 61400) 9A( 25) 96125) 9C(2U) »D1 20) 2AA( 200) 2681200) »CC(200) » 


100(100) ste (i0u) 


GEE = J2e6* EXPE (Oe4*LUGF(G(L))) 

RAR Oo l*EXPF (Oe325*LULE(G(L))) 

V43 RAR#* 2 

IF (V43—-1e0) 41094119410 

RE = (B(19)/B(20) )#(GEE/G(L) )¥#(2e0*RAR/(1¢0~V43) ) 
VA = (B(19)/8H(20))*(GEE/G(L)) 

EE(3) = LeO/SORTFIRE*##2+VAR #2 ) 

EE(4) = ATANFI(VA/RE) 


GO TO 412 
EE(3) = O20 
EE(4) = O00 
RETUKN 


CND 





“tt 
rir) 
' 


jena ae | 





410 


411 
412 


306 
307 


308 


309 
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SUBROUTINE MAI Z(G sA st oC eVolk sAA st »sCCaeDDoeEt ) 
DIMENSION G1400) 0A(25)98B(25) 2C(2Z0) 9D(20) #AA( 200) »9BB1( 200) sCC (200) » 


LDD0( 100) sEE (LOU) 


GEE = (Zeb EXPE (OeG*LOGE (G(L))) 

RAR = Oel#trArl (O0e325#LUGLF ((L))) 

V43 = RAR? 

LF(V43~100) S1Us% Ll oe lu 

RR = (BOELY) 61620) ) (GEES OIL) eC SeO*RARS (1 60-V43)) 

Ree= (CEI 201) ee (OLE/GTE 

GO TO 412 

RR = A(7) 

VA = (B(4) **2#3(5))* (1B 16) 46°00) *#A11))7(816)4814)*A01))) 
= (G(L)*B(3))/7(VA#(B (3) *¥*24+1¢0)) 

VE = RE*8K(3) 


V47 = SUKTFURREHZ+AA#* 2 ) 
V48 = ATANFCAA/RR) 
V49 = SURTr (Re*e24+VE #2 ) 
V50 = ATANF(VE/RE) 


TOP] = V&/#V4Y 

TOP2 = V4stV50 

BOTL]L = SURTF( (KE HKK) ** 2+ (AAFVE ) *#2 ) 
BOT2 = ATANE ((AAHVE)/(RE+RR) ) 


EE(3) = TUPI/BOTL 
EE(4) = TOP2-BOl? 
RETURN 

ENO 


THIS 1S FOUNDATION 

MATRIX THkéet 

SUBROUTINE MAI3 (UsAstsCeDol sAA stb eCC aL eEE ) 

DIMENSION G6140VU) sA025) 96125) 9C120) 90120) 29AA(2Z0U) 9681200) »CC(2UU) » 


1D0(100) sEE (100) sVB( 20) 


V43 = LeO/((7) *(L)) 
DO 405) [eames 

LF (DC1L+10) **2-G(L)**2Z) 30153025301 

VB(I+10) = (GILIFCELI*CCLFLO) I /SIDET 410) #*®2—-G(LJ #*2) 

GO TO 305 

IF (G(LP*CCLI*COCL+10)-0540) 30493039304 

VB(1+10) = 0020 

GO TO 305 

VB(I+10) = AC?) 

CONTINUE 

V78 = VBCLIIFVE LIZ) FV O13) FV0014)4Vb015)4+VB(616)4VB017)+Vb(18)+ 


1VB(19)+vVvB(20) 


V79 = V7B/B(7) 

V22S = VI9-V43 
IF(V229=0020) 30693089307 
EE(6) = A(6) 

GO TO 309 

EE(6) = Al(2) 

GO TO 309 

EE(6) = Oce0 

FE(S) = AC?) 

GO TO 31lé2 

BE(3) = eu /ABSE(Ve?c9) 
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312 RETURN 
END 


THIS IS FOR FOUNDALION Whit DAMPER FOLLOWING 
MATRIX THREt 
SUBROUTINE MAES (GeNottelelolenAoststs all sbboett ) 
DIMENSION G1400) sAl/o) ehs(25)9CF 20) 90020) SAAC ZUU) #bB1 200) CC (200) » 
LOD(100) sEE (LOO) eVEI20) 
V43 = LeO/(BE/)¥G(L)) 
DOe30551 = 1Leld 
PEO +1Ove F2-G(L)* 7) SOND Ua oe 
301 VBC(IF10) = (GILI*FCEIPIYXCtL+ LOIS (001410) FX 2-G(L) ¥¥*2) 
GO TO 30> 
302 IF (G(LIFCOCLIACCIT+10)—-00e0) 304935039304 
303 vBl(It+10) = Q.0 
GO [TO 305 
3046 VB(It10) = AC?) 
305 CONTINUE 
V78 = VBO11)I4VB012)4VH013)4V0014)4VKb015)4+VB010)4+VB(17)4+VK(18)+ 
1VB(19)+VbB(20) 
V79 = V@b/BI(/) 
V229 = VI9~-V43 
IF(V2297"00e0) 30693049401 
306 EE(6) = Al6) 
’ G@ FO 30Y 
30/7 EE(6) = Al(2) 
GO TO 309 
308 EE(6) = Oe«0 
EE(5) = At7) 


GO TO 312 
309 t&05) = LeO/ABSIFIV22Y) 


312 / (16) = LeV/JA( LY) 
EE( 14) = LeOQtlEt (5d) EE lo) ) 
VA 9s HER CS)- Lelio) ) «so lit VE 
RE = LeOF(EE(S) #EbL( 16) )*cCUOSF (ck (6)) 
EE(15) = ATANEC VASE) 
RETURN 
END 


FOUNDATION With DAMRER VakiaAbbe Wilti th EQVUENCY 

MATRIX THREE 

SUBROUTINE MALS (CSOsAdstss: stroll oA srstss CC stUOLL ) 

DIMENSION GL4&00) 9A Zd) Ob 25) 9C 120) eDI 20) sAA( 200) »9BBI 200) »CC( 200) » 

1DD(100) sEE (LOU) eVBI 20) 

V43 = Le0U/( BE) UCL) ) 

DO 305 | = Jel 

TRIOC IT +L Uy RU 2=-G(L IX #2) SUle 407 ea 
201 VB(It10) = (G(L)#FCOL)*C11t10)) fCo01+10)** 2-ol(L) **2) 

GO TO 305 
302 IF (G(LIACCL)FCOCE+L0I~GOe0) 30493039304 
303 VB(It+10) = U-0 

GO TO 305 
304 VB(1+10) = At /) 
305 CONTINUE 

V7B = VtsldldeVtslaZ2dtVaol LAV tVish lad tvBblyodtVE Clo +VBO1L/)+VBEL8) + 
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Lvi(lvb+Vis( 20) 

Viy = Vilu/bl/) 

V229 = VIV-V43 

IF (V229-U0e0) 305930K 9307 
306 EE(6) = Al(6) 

GO TY 30Y 
307 EE(6) = Al2) 

GO TO 309 
308 EE(6) = Oe 

EE(5) = A(7) 

GO TO 312 
309 EE(5) = LleV/AKSFIV229) 
312 GEE = 72e6*e XE (DOe4#*LULE (O(L))) 

RAR = Oel*EXPE (Oe 320* LUGE (U(L))) 

V44 = RARHK#?Z 

IF (V44—-12e0)5002501 9500 
506 RE = BIL)I/SA(1LO) #(GEE/G(L) )#(2e¢O*RAR/(1007V44)) 

EE( 16) = LeO/KRE 

GO TO 502 
501 Et(16) = OeD 
502 EE(14) = 1LeO+(eCbE(D) FEE( 16) ) 

VA = (EE(S) *#EER(16))#SINE(EE(6)) 

RR = LeOF(EE (DS) #**EE(16))*COSE (EEO) ) 

EE(15) = ATANE(VA/RER) 

RETURN 

END 


THIS IS COLUMN HULL MATRIX G 
SUBROUTINE MATS (GeAsbeCeDol epAAstbb a CC eV ett) 
DIMENSION G6(400) sA(25) 9B(25) 9C (20) 2D( 20) eAA( 200) #BB( 200) »CC( 200) » 
1DD(190) s£& (100) »VB( 20) 
V43 = 1e0/(8(8)#U(L)) 
DO 405 [1 = led 
IF (DD( 1) **2-G(L)*#2) 40104022401 
40] VBI) = (GILI#F2e0* DUI IT +5) I / (bOI *(OO01) **¥2—-G(L) ¥*2) ) 
GO TO 405 
402 VB(I) = A(7) 
405 CONTINUE 
V78 = VB(1)+VB(3)4+V8B(5)4+V43 
V79 = VB(2)+VB(4) 
EE(7) = SQRTF(VIB##24+V79%#2 ) 
IF (V79-020)40704060407? 
406 EE(8) = A(6) 
407 €F&(8) = ATANFI(V78/V79) 


V45 = BI9) +(B(11)#GO(L))/A(3)4A(3)*8(12)/G6(L) 

V46 = 36620/6(17) *G(L)I-b(LI/B(16) 

V45 = V4o4+bl lo) 7010046015)" 2: V46%%2) 

V47 = BO10)46(14) *G(LI7 s86 00-019) /G(LI 416015) F¥2KV4G6) (1 eOtH( 15) 


L¥2*V46*2 ) 

EE(9) = SORTFIV4S*#24V4G47# XQ) 
EE(10) = ATANFI(V4//V45) 

V50 = CECT) *EE(9) 

V51 EE(8)+EE(10) 

V52 V>0*COSR(V>l} 

V53 VSO*SINF(V51) 





101 
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V52 = 1leOQtV52 

EE(11) = SQRTF(V52**24+V539%*2) 
EE(12) = ATANF(V53/V52) 
RETURN 

END 


MATRIX FIVve 

THIS IS WATER 

SUBROUTINE MATS (GsAsBoeCoDol sAAsbbeCCsDUecE) 

DIMENSION G6(400) 29A(25)9B(25) 29C(20) 9120) sAA( 200) 9681200) »CC( 200) » 
100(100) »sEE(100) 

EE(13) = 1-0/B(18) 

RETURN 

END 


SUBROUTINE MATS (GeAsBeCaDol sAAsbbeCCeDD x EE ) 

DIMENSION G(400) 9A(25) 98(25) C120) 9D( 20) 29AA( 200) 98B( 200) »9CC(200) » 
100(100) »€E(100) 

EE(13) = G(L)**2/B(1) 

RETURN 

END 


MATRIX MULTIPLICATION 
SUBROUTINE MAMUL (IeJeK0V 0A) 
DIMENSION V(200) 2A(25) 
T=1 

J29 

K=68 

K=Kt2 

V4leV(I)*V(J) 
V422V( 141) 4+V(U+1) 
V43=SINF(V42)*V41 
V442=COSF (V42)*V41 
V45sV(14+2) *V(5U+4) 
V462V(14+3)4+V(U+5) 
V47T=S INF (V46) FV45 
V4B=COSF (V46)#®V45 
V49=V444+V48 
V502=V43+V47 

V(K) =SQRTF(V50**24+V49 #2 ) 
IF (V44+V4B8) 50225950 
V(K+1) = ATANF(V50/V49) 
GO T9 26 

IF(V43+V47) 51s 279 51 
V(K+1) = A(6) 

GO TO 26 

V(K+1) = Oe0 

IF (K-72) 10999 102s 1100 
IF (K-74) 10999 3s 1101 
IF (K-76) 10999491102 

IF (K-78) 10999591103 

IF (K7-80) 10992621104 

IF (K-82) 10999791105 

IF (K-84) 10999891106 

IF (K-86) 10999991107 





Py 
11¢8 
1109 
1110 
Pa 
2112 
Pid3 
1099 


102 


10 


FL 


iZ 
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IF (K-58) 109991021108 

IF (K-90) 10999111109 

IF (K-92) 109991291110 

IF (K~94) 109991341111 

IF (K-96) 109991491112 

IF(K~“98) 109991591113 

IF(K-100) 109991651099 

I = 1+16 

J = Jtl1l6 

GO TO 101 

V102 = (V(70)*#V(72) D¥SINF(VO71)4+V073)) 
V103 = (V(70)*V(72))* COSF(V(71)4+V(073)) 
I = I1~16 

J = J-14 

GO TO 101 

I = 1+20 

J 2 J+14 

GO TO 101 

V104 = (V(74)*V(76))% SINFC(V(75)4V(77)) 
V105 = (V(74)#V(76))*COSF (V7 DtV(77)) 
V(106) = V104+V102 

V(107) = V105+V103 


I] = I~16 

J = J~16 

GO Tl 02 

I = I+1l2 

J = Jt+16 

GO TO 101 

Vl102 = (V(78)*V(80))*SINF(V(79)+V(81)) 
V103 = (V(78)#V(80)) *COSF (V(79)+V(81)) 
I = I-12 

J = J-14 

GO TO 101 

I = [+16 

J = Jt+14 

GO TO 101 


V104 = (V(82)*V(84) ) #SINF (V(83)4+V(85)) 
V105 = (V(82)#*V(84))*COSF (V(83)+V(85)) 
V(108) = V1l04 +V102 
V(109) = V1l05+V103 


I = I-16 


VIO4=(V(90)*V(92) )* STINE IVI9L) FV 095) ) 





13 


14 


i? 


1? 


22 


28 
53 


30 
29 
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V105 = (V(90)*V (92) )*COSF(VI91L) tVI93)) 
V(1l10) = V1l04 + Vive2 
Vii) = Vi0or V1US 


} =21=16 
J = J-18 
GOTO S 04 
I = Itle 
J = J+168 
GO TO 101 


V102 = (V(94)*V(96) D*SINFE (VI95) +V097)) 
V103 = (V(94)#V(46)) KCOSF (VI95)4+V(97)) 


I= i= le 
J = J=-16 
GO TO 101 
I] = 1+16 
J = Jt16 
GO TO 101 


V104 = (V(98)*V(100) )* SINE (V(I99)4V(101)) 
V105 = (V(98)*V(100)) #COSF(V(99)4+V(101)) 
V(ll2) = V104 +V102 
V(113) = V1O5+V103 


I = I1+85 
J= J+6 
K = K+20 
K = K+2 


V43 = SORTF(V(1L) #*24+V C1 +) )¥*2) 
IF (VA I1t1) +000) S2e25ede 

V4G = ATANFIV(L)/VG141)) 

GO TU 29 

LF (V(1)4+00¢0) 53930953 

V44 = A(6) 

GO TO 29 

V44 = 0e0 

V(K) = V43*V(J)*SINE (V444V(I41)) 
V(K+1) = V43#V(IS) #COSFI(V444+V(JI+1)) 
IF (K-124)1098 01891115 

IF (K-126) 109801901116 

IF (K~-128)109892001117 

IF (K-130) 109802191118 

IF (K~132) 109802291119 © 

IF (K-134)109802391120 

IF (K-136) 109892491098 


1 = 1+4 
Jo= J+4 
GO TO 17 


V138 = V(122)+V(124) 
V139 = V(123)+V(125) 


I = [2 
J = J=4 
GO TO 17 
1 = [+4 
Jaewo 44 
GO TO 17 


V140 = V(126)+V(128) 
V141 2 V(1l27)+V(129) 
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l=I-5 

J=J~-2 

GO TO 1? 

T=I+4 

J=Jt+4 

GO TO 17 
V142=V(130)+V(132) 
V143=V(131)+V(133) 
lsI~-2 

J=J~4% 

GO TO 17 

l=I1+4 

J=ZJ+4 

GO TO 17 
V144=V(134)+V(136) 
V145=V(135)4+V(137) 
RETURN 

END 


THIS IS TRANSLATION 
SUBROUTINE TRANS (V.sEE) 
DIMENSION V(200) 0EE(100) 
V(1l) = 1leA 


V(2) = Oe0 
V(3) = EE(1) 
V(4) = EE(2) 
V(5) © Oe 
V(6) =0e0 
V(7) = 1leO 
V(8) = Oed 
V(9) = 1¢0 
V(10) = 020 
Vill) @ OeO 
V(12) = Oed 
V(13) ='EE(3) 
V(14) = EE(4) 
V(15) = 1¢0 
V(16) = OeO 
V(17) = 140 
V(18) = OeO 
V(19) = EE(5) 
V(20) = EE(6) 
V(21) & Oe 
V(22) = OeO 
V(23) 2 1¢0 
V(24) = Oed 
V(25) = 1¢0 
V(26) = 040 
V(27) = EE(9) 
V(28) = EE(10) 
V(29) = EE(7) 
V(30) = EE(8) 
V(31) = EE(11) ‘ 
V(32) = EE(12) 
V(33) 2= 120 
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V(34) = 020 
V(35) = 040 
V(36) = 020 
V(37)) = EE(13) 
V(38) = 0Oe0 
V(39) = 1¢0 
V(40) = 020 
RETURN 

END 

FORTRAN 


THIS IS TRANSLATION 
THIE IS FOR FOUNDATION WITH DAMPER FOLLOWING 
SUBROUTINE TRANS (Voekt) 


DIMENSION V(200) eEE( 100) 


Vt(l) = 1 
V(2) 
v(3) 
V(4) 
V(5) 
V(6) 
V(7) 
V(8) 
V(9) 3 
v(10) 
Vill) 
vt(l2) 
V(13) 
V(14) 
V(15) 
V(16) 
V(i7) 
Vv(18) 
V(19) 
V(20) 
v2.) 
v(22) 
v(23) 
V(24) 
v(25) 
V(26) 
Vi27) 
V(28) 
Vv(29) 
V:030) 
V(31) 
Viz 
V(33) 
V(34) 
V(35) 
V(36) 
V(37) 
V(38) 
V(39) 
V(40) 


keuman @ 
>) 


out ok bt nt 8 We Hh Rhee ket one OO oka ak ae Oa ae 


00 


020 
120 
0.0 
1.20 
020 
EE(9) 
EE(1TO) 
BE (CT? } 
EE(8) 
EE(11) 
EE(12) 
120 
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RETURN 
END 


SUBROUTINE ANSW (AsG1lsGoE ol oV) 

DIMENSION A(25)9G61(400) 9G(400) 9E (400) »V(200) 
V138 = V(1l22)+V(124) 

V139 = V(123)+V(125) 

E(L) = (A(5)/A(4)) *#LOGF(1¢O0/SQRTF(V1I13B8**24+VI99*K2 ) ) 

Gl(L) = G(L)I/A(3) 

G2 = LOGF(G1(L))/A(4) 

MFREQ = (((G2—-120)*92000)/3e0)+10200 

MT = (((E(L)+7000)*87000)/14000) 415200 

CALL BRIPIC(MFREQoMT ) 

CALL BRIPIC (MFREQ MT) 

RETURN 

END 


SUBROUTINE CHOMEG(E sL9G9A) 
DIMENSION E(400) 9G(400) sA(25) 
UeCGULY=AL9) ) 251925900200 

V502 = ABSF(E(L)~-E(L~-1)) 

V1012 = ABSF(LOGF(G(L))-LOGF(G(L=-1))) 
1F(V502=540)2549253 9223 

IF (V502-225)25092559255 

IF (V1012-0202878) 25092539250 
G(L+1) = EXPF(LOGF(G(L))+0e02878) 
GO TO 252 

G(L+1) = EXPF(LOGF(G(L))+-e05757) 
RETURN 

END 


SUBROUTINE PRINO (GleGsEosL) 

DIMENSION G1(400) »G(400) sE (400) 

PRINT 2209 (G(I)sGl( IT) sE(T)sT2l ol) 

FORMAT(8H OMEG = F8els9H FREQ = Flels6H T = E1306) 
CALL FRAME 

RETURN 

END 





(89 ) 


APPENDIX B 
PRESENTATION OF RESULTS 


The following photo plots were produced by the Cathode Ray Tube 
Recorder of the IBM 709 Digital Computer. The computed points were 
connected by solid lines by the authors to clarify the plots. 


Due to the space limitations of the scope presentation the 
abscissa or frequency scale is scaled in integers for three decades. 
The first decade on all of the plots extends from 10 to 100 cps, the 
pecond decade from [00 to 1000 cps) asudmthesthird trom 1000 to We 7es> 
eps. The ordinate presents transmission in decibels varying from 
-70 to +60. 
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